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1 Introduction 

At low temperature T, a significant difference between the behavior of crystals 
on the one hand and disordered solids on the other is seen: sufficiently strong 
disorder can give rise to a transition of the transport properties from con- 
ducting behavior with finite resistance R to insulating behavior with R = oo 
as T — > as was pointed out by Anderson in 1958 pi. This phenomenon is 
called the disorder-driven metal- insulator transition (MIT) an d it is 

characteristic to non-crystalline solids. The mechanism underlying this MIT 
was attributed by Anderson not to be due to a finite gap in the energy spec- 
trum which is responsible for an MIT in band gap or Mott insulators S. 
Rather, he argued that the disorder will lead to interference of the electronic 
wave function ip(x) with itself such that it is no longer extended over the 
whole solid but is instead confined to a small part of the solid. This localiza- 
tion effect excludes the possibility of diffusion at T = and thus the system 
is an insulator. 

A highly successful theoretical approach to this disorder-induced MIT 
was put forward in 1979 by Abrahams et al. ra. This "scaling hypothesis 
of localization" details the existence of an MIT for non-interacting electrons 
in three-dimensional (3D) disordered systems at zero magnetic field B and 
in the absence of spin-orbit coupling. The starting point for the approach 
is the realization that the sample-size (L d ) dependence of the (extensive) 
conductance 



should be investigated |3| with a denoting the conductivity, d the spatial di- 
mension and g the dimensionless conductance. On the other hand, for strong 
disorder, the wave functions will be exponentially localized with localization 
length A and thus the conductance in a finite system will be 



G = o-L 




(1) 



g ~ exp(-£/A). 



(2) 



Defining the logarithmic derivative 



ding 



(3) 



dlnV 
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we see from Eq. (0) that < for d = 1 and 2 and thus an increase in L 
will drive the system towards the insulator, there are no extended states and 
no MIT. However, the (3 curve for d = 3 is positive for large g and negative 
for small g and an increase of L will drive the system either to metallic or to 
insulating behavior. 

The scenario proposed by the scaling hypothesis is that of a continuous 
second-order quantum phase transition Q] . Then in the vicinity of the critical 
energy E c the DC conductivity a and the localization length A should behave 
as 

a oc (E ~ E c ) s for E > E c , (4) 
A cx (E c - E)~ v for E < E c . (5) 

with s — u(d — 2) due to further scaling relations ||. Introducing a similarly 
defined dynamical exponent z for the temperature scaling as o~(T) oc T 1 / 2 , 
we can write the full finite-temperature scaling form as 



a( M ,T)cx [(p-E c )/E c ] s f 



T 



[(fj. - E c )/E c 



(G) 



with \i the chemical potential, T the scaling function and z = d The 
special energy E c is called the mobility edge Q and separates localized states 
with > E c from extended states with \E\ < E c . States directly at the 
transition with E = E c are called critical and will be examined later in much 
detail. In this way the disorder-driven MIT has been reformulated in terms 
of the theory of critical phenomena 



2 Experimental evidence in favor of scaling 

Much work has subsequently supported these scaling arguments at B = 
experimentally, analytically and numerically ||. The MIT can be observed 
by measuring the conductivity on the metallic side and the dielectric sus- 
ceptibility on the insulating side of the transition B. For doped Si:P, many 
experiments have been performed following the original work of Paalanen 
and Thomas ||. The one-parameter scaling hypothesis has been beautifully 
validated in these experiments by, e.g., constructing scaling curves for the 
conductivity |1Q]. The recent experiments in Si:P 0,12, 13, l(|are concerned 



with the exact estimation of the critical exponent v as in Eq. (W). The current 
estimate is v « 1. The interest in the exact value of v arises since compensated 
semiconductors apparently have u « 1 as do amorphous metals 14,1^, p|. On 



the other hand, for uncompensated semiconductors one had previously found 
v k. 0.5 &p|. The recent estimations of v w 1 for uncompensated Si:P pfl |, 
based on a careful scaling analysis according to Eq. (^|) and a considera- 
tion of various temperature regimes, may suggest at last a resolution of this 
"exponent puzzle" pq. 
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Other experiments in 3D have been performed, e.g., on Si: B Scal- 
ing according to Eq. (||) yields v — 1.6. The large value of v — as compared 
to the Si:P data — was attributed to the presence of interaction effects. An 
experimentally convenient way to construct very homogeneously disordered 
samples is the transmutation doping technique which uses the ho- 

mogeneous properties of neutron rays. Recent scaling results then suggest 
v = 1.6 ±0.2 Q. 

As the localization phenomenon in disordered solids is intrinsically due 
to the wave nature of the electrons, it can also be observed in other systems 
exhibiting wave motion ||] . Localization has been studied, e.g. for water waves 
[ p3| in shallow basins with random obstacles, for light waves |24|, ^5[]2(| in 
the presence of a fine dust of semiconductor material, for microwaves p7| . [28| 
in microwave cavities with random scatterers, and also for surface plasmon 
polariton waves |p9[ on rough semiconducting surfaces. 



3 Scaling and the Anderson model of localization 

In order to describe a disordered system, let us consider the Anderson model 
of localization [0 , 

\ja)(m+J2^\ja)(ja\. (7) 

The off-diagonal matrix tj a ^p denotes the hopping integrals between the 
states {a} at sites {j} with the states {(3} at sites {k} and represents the 
discretization of the kinetic energy. For simplicity, one usually assumes that 
a = P = 1 such that there is only one state per site. Moreover, the hopping 
is usually restricted to nearest-neighbor sites. The disorder is incorporated 
into the diagonal matrix tj a , whose elements are random numbers usually 
taken from a uniform distribution [— Wy2, W/2] with W parameterizing the 
strength of the disorder. Other distributions such as Gaussian and Lorentzian 
have also been investigated ]30|,|3^,|3l], ||l [33"t . 

This model has been used extensively in conjunction with powerful nu- 
merical methods in order to study the localization problem ||. MacKinnon 
and Kramer |3^,|36| have numerically verified the scaling hypothesis by show- 
ing that one can find the scaling behavior outlined in section [I[ For 3D the 
corresponding scaling curves have two branches corresponding to the metal- 
lic and insulating phases, whereas in ID and 2D only the insulating branch 
exists. Recent numerical results — some of which shall be presented in the 



coming sections — indicate that v — 1.58 ± 0.02 |37, 3§|] which is in excellent 



agreement with the newer experiments reviewed above. 
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4 Numerical methods and finite-size scaling for 
disordered systems 

The preferred numerical method for accurately computing localization lengths 
in disordered quantum systems is the transfer-matrix method (TMM) ]39|,|35|, 
E36LEJ. The TMM is based on a recursive reformulation of the Schrodinger 
equation such that, e.g. in a 2D strip of width M, length N 3> M and uniform 
hopping tj t k = 1 between nearest neighbors only, 

1pn-\-l,m — {E £n,m)lftn,7n ^n,m — 1 V^n.m+l lfin—l,m ($) 

where ip n ,m is the wave function at site (n,m). Eq. (^) can be reformulated 
into a matrix equation as 

(ip n+ i\ (E-e n -H ± -1 
\ i>n ) \ 1 

where ip n = (4>n,ii ■ ■ ■ >ipn,M) T denotes the wave function at all sites of the 
nth slice, e n = diag(e„ > i, . . . , e n ,M), an d H± represents the hopping terms in 
the transverse direction. The evolution of the wave function is given by the 
product of the transfer matrices tjv = T/vT/v-i • • ■ T 2 Ti. Strong fluctuations, 
which increase exponentially with the system size, govern the evolution of the 
wave function and thus the behavior of the transmission coefficient through 
the samp le |46|^2|,|4^,^3|,|4l| , ^4| . Only the logarithm of the transmission coeffi- 
cient |||,^|^,^0|,|42j,^lJ and the logarithm of the conductance [|6[|^,|8| are 
statistically well-behaved self-averaging quantities. According to Oseledec's 
theorem [|9| the eigenvalues exp[±7i(M)] of _T = limjv^ oc (rj vr TAr) 1 / 2JV exist 
and the smallest Lyapunov exponent 7 m i n > determines the localization 
length X(M) = 1 / 7 m i n at energy E. The accuracy of the A's is determined as 
outlined in Refs. ^,^6) from the variance of the changes of the exponents 
in the course of the iteration. Usually the method is performed with a com- 
plete and orthonormal set of initial vectors (i[>i,ipo) T . In order to preserve 
this orthogonality, the iterated vectors have to be reorthogonalized during 
the iteration process. 

For small disorders, the A(M) values are of the same order of magnitude 
as the strip width M and thus subject to finite-size modifications. In order 
to avoid simple extrapolation schemes, a finite-size scaling (FSS) technique 
had been developed in Refs. |55[[36|,^(| based on real space renormalization 
arguments for systems with finite size M and intimately related to the orig- 
inal scaling approach |^,|5^]. The connection to the experimentally perhaps 
more relevant finite-temperature scaling as in Eq. (|J) is based on the idea 
that a finite system size M may be assumed to be equivalent to a measure- 
ment at finite temperature T since a finite temperature induces an effective 
length scale beyond which the electrons will scatter inelastically and thus lose 
the phase coherence necessary for quantum interference |^,pq|. Scaling the 
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\(M)/M data for various values of W onto a common scaling curve, i.e., 

A(M)/M = f((/M). (10) 

is the analogue of Eq. (||). One determines the FSS function / and the values 
of the scaling parameter £ by a non- linear Levenberg-Marquardt fit |pq| , 
see also T. Ohtsuki's contribution in this volume. For diagonal disorder in 
3D, this scaling hypothesis of localization has been shown to be valid with 
very high accuracy, and the £ values of the extended (localized) branch are 
equal to the correlation (localization) length in the infinite system. A similar 
method based on the recursive Green's function technique is discussed by A. 
MacKinnon in this volume. 

A number of more indirect numerical approaches to the MIT have been 
developed that either only require selected parts of the spectrum — so-called 
energy-level statistics (ELS) — or a few selected eigenvectors in the spectrum 
- so-called wave- function statistics (WFS) [|5l|,|54 j53||. These methods are 
based on the connection of Anderson localization to random matrix theory 
(RMT) @,||. 

The (inverse) participation number |56| represents another possibility to 
distinguish localized states (57[, its scaling with the system size yields a char- 
acteristic fractal dimension 64 58|,[59|, |63| , |60| , |64| , |6l| | . A generalization to higher 
moments of the spatial distribution of the wave function leads to the multi- 
fractal analysis (MFA) |6^,^], where one computes a spectrum of exponents 
to describe scaling properties of the wave function ]6^,[7^,|7l|,|7^,|7r^,[7^, [75| , |6^ , 
[38], |6{], f70|. For a given disorder, one can then read off from the system-size 
dependence whether the spectrum tends towards the metallic, insulating or 
truly critical behavior. At the MIT the singularity spectrum of the MFA is in- 



dependent of the system size 76 , 75|,|77|] . Thus we again have a means of study- 
ing the MIT. Note that WFS and MFA are not independent of each other, 
but one may be derived from the other j78) . The fractal characteristics of the 
eigenstates, i.e their scale incariance, can be beautifully visualized (7|,|8(| by 
displaying the curdling of the wave functions at the MIT. 



5 Scaling for non-interacting, disordered systems 

5.1 Early scaling results for the isotropic Anderson model 

Following the seminal papers of MacKinnon and Kramer j3f], the TMM 
and FSS approach was used to determine the phase diagram of localization 
i.e. the MIT in the entire (E, W) plane for the isotropic Anderson model with 
diagonal box, Gaussian and Lorentzian disorders |34|,|31, 32, 33L ^lj as well as 



at the band center for a binary |82J] and a triangular distribution |83j]. In 
2D, not only the usual square lattice, but also a triangular and honeycomb 
lattice was studied Q]. The derivation of the critical exponents in these 
early studies was, however, impeded by the relatively large error of 1 % of 
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the raw data and the limited cross section size of typically up to M 2 = 13 2 
Much effort was needed to determine reliable values of 




We note that FSS has also been successfully employed for the analysis 
of ELS data and the derivation of critical exponents from the cumulative 
level-spacing distribution and Dyson-Metha statistics p7|,|8S|]. 

The scaling of the participation number has further been investigated 
for the Anderson model on 2D and 3D quasiperiodic lattices, too H,§9[. By 
means of FSS of the participation number, the MIT and the critical exponent 
could be computed [jS9| . 

The entire phase diagram of localization has also been established by MFA 
of the isotropic Anderson model with box, Gaussian and binary disorder 
[[75), determining the parameter combinations of E and W for which the 
singularity spectrum is scale-invariant. 



5.2 The Anderson model with anisotropic hopping 

As shown in section [j], there is no MIT in 2D in the absence of many-body 
interactions, magnetic field and spin-orbit interactions. Furthermore, the 2 + e 
expansion within the non-linear a model Q] and numerical studies based on 
TMM data for bi-fractals (9lJ suggest that the critical exponent in 2 + e 
dimensions changes continuously as e — > for e between and 1. Thus one 
is led to ask the question whether a similarly continuous change does also 
happen, if we vary the hopping elements anisotropically (see C. Soukoulis 
in this volume). E.g., we decrease the hopping homogeneously in one or two 
directions yielding weakly coupled planes or chains. This then might model 
a transition from 3D to 2D or ID, respectively. 

We have investigated this problem using MFA ||,||, TMM H,|| and 
ELS |)4|,|95|,|88| together with FSS. We find that the critical disorder changes 
continuously with decreasing hopping strength t a such that W c ~ i", where 
a is close to j for planes and ^ for chains. Here a represents z for coupled 
planes and y and z for chains. In Fig. |] we show examples of FSS in this 
model. Note that the small relative error 0.07% of the raw data as well as 
the large system cross sections up to M 2 — 46 2 |)4|] (not shown here) for 
TMM make it possible — besides taking into account non-linear deviations 
from 1 1 — 1 — to determine estimates for irrelevant scaling exponents | Bq | . 
The value of the critical exponent v is not affected by the anisotropy |)6| and 
retains its usual value 1.6 ± 0.1 as in 3D [^7[^8|. Thus the 2D and ID cases 
are reached only for t a — and we observe the 3D MIT at any finite t a - 

The ELS and the singularity spectrum of the MFA at the MIT are inde- 
pendent of the system size and this size independence can be used to identify 

the MIT QgHSHSil. However ' both ELS and MFA are influenced 
by the anisotropy and change considerably in comparison to the isotropic 
case. Also, the eigenfunctions are different from the isotropic case. Therefore 
ELS and the MFA singularity spectrum at the MIT are not universal, i.e., 
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Fig. 1. Left: Reduced localization length for coupled planes with t z = 0.1 and 
M = 5, 6, . . . , 17. The lines are fits of the data according to |3$| with n r = 3 and 
m r = 2. In the inset we enlarge the central region without the data points to 
show the shift of the crossing point. Right: Scaling function and scaling parameter, 
shown in the inset, corresponding to the fit in the left panel. The symbols distinguish 
different W values of the scaled data points. 



not independent of microscopic details of the system. This dependence on 
microscopic details is similar to the dependence on boundary conditions es- 
tablished recently for the ELS at the MIT JT?5o| , [lOl] ] . However, the critical 
exponent should, of course, be universal. As it turned out [|4[^], FSS can 
be applied successfully to various statistics of the spectrum, most accurately 
to the integrated E% and A3 statistics. In Fig. g, we give scaling results for 
the integrated £2 statistics (77) and its derivative, the so-called level com- 
pressibility Xi both of which have been computed from spectral data with 
error 0.2 - 0.4%. The critical exponents v n = 1.43 ± 0.13, v x = 1.47 ± 0.10 
derived from these data p5[ |, although less accurate, are in agreement with 
the above-mentioned values. 



5.3 The Anderson model with random hopping 

Let us change the model (fy such that all nearest-neighbor t values can be 
chosen randomly [102] with, e.g., t 6 [c — w/2, c + w/2] with c and w denoting 
the center and width of the distribution. In this parameterization, the ordered 
tight-binding model is recovered in the limit c — » 00 after a suitable rescaling. 
The DOS has a peak at energy E = for any strength of hopping disorder 
— known in ID as Dyson singularity [ 103 1 — and the localization length at 
E = diverges even in ID [ 1 04 1 . 



We have studied the random hopping model in 2D p^, 10q , [l07[ by TMM 
and by direct diagonalization of the Hamiltonian matrix and have shown 
that the singularity in the DOS still exists for bipartite square lattices where 
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Fig. 2. Left: The one-parameter scaling dependence of r/ on £ for different system 
sizes N and disorders W £ [6, 12]. Right: The one-parameter scaling dependence of 
X- The dashed line indicates the value Xc = 0.27 at the MIT obtained from this fit. 



the energy spectrum is strictly symmetric around E — 0. Furthermore, the 
localization length is also diverging at E = |108 . 109,113, 11 1 1 , see also S. 
Evangelou's contribution in the volume. 

For sufficiently large energies, it has been suggested [11 1 that the diver- 
gence of the localization length at the band center may be described by a 
power law, whereas it takes more complicated form below a certain crossover 
energy E* . Our results at 0.1% error in the raw data for 2D suggest that the 
localization lengths exhibit power-law behaviour in a wide energy range with 
lower bound E m i n « 10~ 7 and non-universal exponents v ss 0.25 [112], see 
Fig. |^. For smaller energies we observe some deviations, however, there is also 
a possibility that this may be an effect of pronounced convergence problems 
which appear for strong hopping disorder. 

In 3D, the hopping disorder alone is no longer sufficient to localize all 
states 1 102] as happens for (uniform) diagonal disorder at W c — 16.5. Results 
of FSS for the 3D system at 0.1% error indicate that the critical exponent v 
is the same regardless whether we study the MIT as a function of E or as 
a function of additional diagonal disorder W. Taking into account irrelevant 
scaling terms, we find that v = 1.59 ± 0.05. Thus the result s are aga in in 
agreement with the usual 3D case and the scaling hypothesis [114. 



5.4 Thermoelectric transport coefficients in the Anderson model 

The conductivity a is the quantity which is most often studied in transport 
measurements of disordered systems [|l7|, [l^, 0, 0, p2L |l3 , 115 , 1(| . However, 
other transport properties such as the thermopower S 1 , the thermal conduc- 
tivity K and the Lorenz number Lq have also been measured [ 116| , 117, lit] . 
In Refs. 11£, 12C, 121], we have studied the behavior of S, K and Lq by 
straightforwardly calculating the integrals in the linear response formulation 
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Fig. 3. Scaling parameter £ vs. energy. Left panel: square lattice, Gaussian t dis- 
tribution, right panel: honeycomb lattice, logarithmic t distribution; filled symbols: 
results for M = 50-100, open symbols: results for M = 110-160 (Gaussian distri- 
bution) or M — 100-150 (logarithmic distribution). 
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Fig. 4. Scaling of thermoelectric transport properties where t = |1 — Ef/E c \. The 
different symbols denote the relative positions of various values of the Fermi energy 
Ef with respect to the mobility edge E c . 



of Chester-Greenwood-Kubo-Thellung |12|,|l2§|2||. The only additional in- 
gredients in our study were an averaged DOS and (a) the assumption of <j(E) 
as in Eq. (Q) j3(], |7j| or (b) an energy-dependent conductivity a obtained 
experi mentally . For (a), we can show that the previous analytical considera- 
tions Jl25j|l26| , |l27|l apply in limited regimes of validity. Thus S/T diverges 
as T — > when the MIT is approached from the metallic side, but S itself re- 
mains constant at the MIT. For (b), we show that the temperature-dependent 
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a, the thermoelectric power S, the thermal conductivity K and the Lorenz 
number Lq obey scaling as shown in Fig. |J 



6 Scaling for interacting, disordered systems 

The research presented in the last section clearly supports the scaling hy- 
pothesis of localization for non-interacting electrons. However, real electrons 



of course interact [128], and their interaction is of relevance for the transport 
properties of disordered systems |129 , 13C , 131], especially in 2D and ID [ 132 1 
where screening || is much less efficient than in 3D. Recently, these theoreti- 
cal considerations received a lot of renewed attention due to the experimental 
discovery of the 2D MIT (jj§ [l3l| . 

In order to theoretically study the effects of the interplay between disorder 
and interactions, in principle one has to solve a problem with an exponentially 
growing number of states in the Hilbert space with increasing system size. 
At present, this can be achieved only for a few particles in ID and very few 
particles in 2D [134,135,136]. However, in 1994 Shepelyansky [137| proposed 
to simply look at two interacting particles (TIP) in a random environment. 
He showed that two particles in ID would form pair states even for repulsive 
interactions such that the TIP pairs would have a larger localization length 



A 2 oc U 2 Xi 



(11) 



than the two seperate single part icles (SP) leading to an enhanced possibility 
of transport through the system 138 1. In Eq. ([ll]) the pair energy is E = 0, U 
represents the onsite interaction strength and Ai is the SP localization length. 
Subsequent works have established that an enhancement due to intera ction 
indeed exists, although the details are somewhat different from Eq. ( [Tlj) [139]. 



6.1 Using decimation to study TIP in random environments 

The obvio us failure of the T MM approach to the TIP problem in a random 
potential | 14Cj , 141 , 142 , 143 ] has led us to search for and apply another well 



tested method of computing localization lengths for disordered system: the 
decimation method [144]. We computed the TIP Green function in ID at 
selected energies for 26 disorders between 0.5 and 9, 24 system sizes between 
51 and 251, as well as 11 interaction strengths U — 0, 0.1, . . . , 1. For each such 
set of parameters, we averaged over at least 1000 samples. Furthermore, we 
constructed FSS curves (see Fig. ^) and from these curves computed scaling 
parameters which are the infinite-sample-size estimates ^(U) [ 145 1 of the 
localization lengths. We found ]146| , |147[ that onsite interaction in ID indeed 
leads to a TIP localization length which is larger than the SP localization 
length at E — 0. However, the functional dependence is not simply given 
by Eq. (pd]). Our data follow ^{U) ~ ^(O)' 3 with an exponent f3 which 
increases with increasing \U\ at E — 0. The best fit was obtained when the 
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enhancement £,2(U) / ^{0) depends on an exponent (3 which is a function of 
U |148| . For values of U > 1.5 we observe that the enhancement decreases 
again; the position of the maximum depends upon W reflecting the expected 



duality between the behavior at small and large U [149 1 
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Fig. 5. Left: Finite-size scaling plot of the reduced TIP localisation lengths X2/M 
for U = (O), U = 0.2 (O) and 17 = 1 (□). The data for U = 0.2 (U = 1) have 
been divided by 2 (4) for clarity. Data corresponding to W = 1 are indicated by 
filled symbols. The two curves at the bottom show the same data for U = 0.2 and 
1 and W < 2.5, shifted downward by one order of magnitude for clarity, but here 
the data for W < 1 are scaled with scaling parameters obtained from the power- 
law fits in the right panel. Right: TIP localisation lengths £2 after FSS for U — 
(O). U = 0-2 (O) and U = 1 (□). The solid line represents ID TMM data for SP 
localisation lengths Ai/2, the dashed lines indicate power-law fits. Inset: Exponent 
a obtained by the fit of £2 tx W~ 2a to the data for (7 = 0, 0.1, . . . , 1. 



6.2 The TIP effect in a 2D random environment 



In Ref. 15C] we have employed the decimation method for TIP in quasi-lD 
strips of fixed length L and small cross-section M < L at E — 0. Analytical 



considerations for 2D [138] predict that the enhancement of A2 at E = 
should be 

~u 2 \l 



A2 oc Ai exp 



t: 2 



(12) 



with the SP localization length Ai cx exp(t 2 /W 2 ) in 2D [§§. We found [jl50 | 
that the enhancement is even stronger and as shown in Fig. ^| the scaling 
curves for U > 0.5 have two branches indicating a transition of TIP states 
from localized to delocalized behavior [151]. The scaling curves for U < 0.2 
show a single branch corresponding to localized behavior. The quality of 
the scaling curves is not as good as in the ID TIP analysis [146], due to the 
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Fig. 6. Left: FSS scaling curves (lines) and reduced localization lengths X2/M for 
TIP in 2D as a function of the scaling parameter £2 for U = (o), 1 (□), and 2 
(V). Ri ght: Critical exponent v. The data point (□) for U = 1 represents the result 
ofRef. pi}. 



smaller samples and smaller number of configurations. Nevertheless, our data 
for 51 interaction strengths and 36 disorder values allow us to map the (U, 
W) phase diagram of the TIP delocalization-localization transition and we 
can study how the critical exponent of the localization length changes with 
U in Fig. H We find that for all U G (0,2], the exponent is systematically 
larger than the critical exponent of the usual Anderson transition for non- 
interacting electrons in 3D. Let us emphasize that this transition is not a 
metal-insulator transition in the standard sense since only the TIP states 
show the derealization transition. The majority of non-paired states remains 
localized. 



6.3 The TIP effect close to an MIT 

The numerical effort to study the influence of interaction directly at the MIT 
in 3D is currently prohibitive. This is true even for TIP, since the problem 
is equivalent to a six-dimensional SP system with correlated disorder. For- 
tunately, there is a ID model which exhibits an MIT driven by increasing a 
local potential. This model is known as the Aubry-Andre model ]152| | which 
we extend by an interaction term, i.e., 

M M M 

H = ( C I-,ri+l C o-,n + ^- C + Mti4,n C ff,n + U nm C n ^ C^C^ C m | , 

o",n— 1 er,7i— 1 n,m— 1 

(13) 

where \x n = 2/icos(cm + (3) with a/2ir an irrational number chosen as the 
inverse of the golden mean a/2-K = (Vo— l)/2. (3 is an arbitrary phase shift. 
The cj. n and c ayn are the creation and annihilation operators for an electron 
at site n with spin a ="[",!■ U nm denotes the interaction between particles: 
U nm = U5 nm for Hubbard onsite interaction or U nm = U/(\n — m\ + 1) for 
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long-range interaction. For /x < 1, all SP states in the model with U = 
have been proven rigorously to be extended, whereas for /j, > 1 all SP states 
are localized p^ , |l5^p^ , p!5^ , p!5^ , p^| ,p^|. Directly at ft = fi c = 1 the SP 
states are critical. Thus the MIT is similar to the MIT in the 3D Anderson 
model, but there are no mobility edges. 

The model has been previously considered at U > from the TIP point 
of view in Refs. [159, 160| , 161]. It has been shown that on the localized side, 
the TIP effects persists, i.e., the TIP localization length A2 > Ai. On the 
extended side, it was ar gued that the interaction leads to a localization of 
the TIP states. In Refs. 162 , 163 ], we have used the TMM and decimation 
together with FSS to study the problem. We use up to M — 377 sites for 
the decimation method with at least 1000 samples. Let us emphasize that 
contrary to the problem with TMM for the TIP situation in finite systems, 
together with FSS the TMM approach can be used to give meaningful re- 
sults. However, the computed localization lengths are no longer directly the 
localization lengths of a TIP pair, but rather measure the influence of the 
presence of the second particle on the transport properties of the first. In 
addition to investigating the onsite interacting case, we have also studied 
long-range interactions in Ref. [162]. We find that whereas onsite interaction 
does not shift the MIT from fj, c = 1, long-range interaction might change the 
MIT towards smaller values fi c m 0.92. 

For the quasiperiodic many-body system with nearest-neighbor interac- 
tion at finite particle density [164 , 165 we have recovered the transition at 
He = 1 independent of interaction strength, provided we consider densities 
like p = 1/2 which are incommensurate compared to the wave vector of the 
quasiperiodic potential — an irrational multiple of tt. Thus, the low-density 
TIP case is comparable to finite but incommensurate densities. On the other 
hand, for commensurate densities, we find that the system can be completely 
localized even for /i <C 1, due to a Peierls resonance between the degrees of 
freedom of the electronic system and the quasiperiodic potential. Whereas for 
repulsive interactions the ground state remains localized, we find a region of 
extended states for attractive interaction due to the interplay between inter- 
action and quasiperiodic potential. Thus, the physics of the model at finite 
densities is dominated by whether the density is commensurate or incom- 
mensurate and only in the latter case by interaction effects. 



7 Conclusions 

In the preceding sections, we have presented results of transport studies in 
disordered systems, ranging from modifications of the standard Anderson 
model of localization to effects of a two-body interaction. Of paramount im- 
portance in these studies was always the highest possible accuracy of the raw 
data combined with the careful subsequent application of the FSS technique. 
In fact, it is this scaling method that has allowed numerical studies to move 
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beyond simple extrapolations and reliably construct estimates of quantities 
as if one were studying an infinite system. Of course, this statement is only 
a short and perhaps too short summary of the seminal paper by MacKin- 
non and Kramer (3f|, in which the FSS technique was first applied to the 
localization problem. 
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